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Signature of a Z2 vortex in the dynamical correlations of the triangular-lattice 

Heisenberg antiferromagnet 
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Dynamical properties of the classical Heisenberg antiferromagnet on the triangular lattice 
are investigated by means of a spin-dynamics simulation and an analytical calculation. While 
the model was suggested to exhibit a topological transition driven by the Z2-vortex binding- 
unbinding, weakness of the associated thermodynamic singularity has made it difficult to ob- 
serve the evidence of the Z2 vortex experimentally so far, only some indirect support. Here, we 
show that the signature of the Z2-voTtex excitation and the vortex-induced topological transi- 
tion can be captured from the dynamical spin correlations. In particular, the dynamical spin 
structure factor exhibits a characteristic central peak around the Z2-vortex transition, and this 
central peak will be a fingerprint of the Z2-vortex excitation. 
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1. Introduction 

Recently, geometrically frustrated magnets have at- 
tracted much interest because of their unconventional 
ordering behaviors, including the possible quantum spin- 
liquid state^^ and the novel ordered states like chiral 
ordered state, spin nematic or multipolar state,**' 
spin-gel state^^ etc. Such novel ordered states as well as 
the associated phase transitions are often borne by novel 
excitations inherent to frustrated systems. One typical 
example might be a Z2 vortex, stabilized in a class of two- 
dimensional (2D) Heisenberg magnet with the locally 
noncoUinear spin order, e.g., the frustrated Heisenberg 
antiferromagnet (AF) on the triangular lattice (Fig.l). 

Indeed, the AF Heisenberg model on the 2D trian- 
gular lattice is a typical example of geometrically frus- 
trated magnets and has been extensively studied as 
the standard model of frustrated systems. In recent 
experiments on the triangular-lattice Heisenberg AFs 
S^3/2 NaCr02,^-") S^l NiGa2S4*^"^**^ and S'=l/2 
organic compounds, K-(BEDT-TTF)2Cu2(CN)3i9-2i) ^j. 



EtMe3Sb[Pd(dmit)2 



22,23) 



a spin-liquid like behavior 



without the standard magnetic long-range order was ob- 
served, while all of these compounds exhibit a weak but 
clear transition-like anomaly at a finite temperature. As a 
possible explanation of the observed experimental behav- 
ior, ref. 6 suggested a Z2 vortex as a key ingredient. Ex- 
perimentally observed weak anomaly was then ascribed 
to the Z2-vortex driven topological transition, while the 
observed low-temperature state to the "spin-gel" state 
where the spin correlation length is kept finite while the 
ergodicity is broken topologically.^'^*-' 

If one considers the simplest case of the classical AF 
Heisenberg model with the nearest-neighbor AF cou- 
phng, the ground state is the so-called 120-degrees struc- 
ture. It was demonstrated by Kawamura and Miyashita 
that the model sustained a topologically stable point 
defect characterized by a two-valued topological quan- 




*Pi_mnil arlrlrnQQ- r\\^iihin(i 



cniTi 



Fig. 1. (Color online) Illustration of a Z2 vortex formed by the 
local 120-degrees structure realized in the antiferromagnetic 
Heisenberg model on the triangular lattice. 



ttmr number, a Z2 vortex^^ (Fig.l). They suggested that, 
though the AF spin correlation length remains finite at 
any finite temperature T > 0, a thermodynamic phase 
transition driven by the binding-unbinding of the Z2 vor- 
tices occurred at a finite temperature. 

The remarkable feature of the Z2-vortex transition 
might be the decoupling of the two length scales, a spin- 
wave (SW) correlation length ^sw and a vortex corre- 
lation length ^y, the latter corresponding to the mean 
separation of free Z2 vortices.®'' At the Z2-vortex transi- 
tion temperature T = Ty, the vortex correlation length 
diverges, while the spinwave correlation length £^sw 
is kept finite (Fig. 2). In the high-temperature regime 
T ^ Ty where ^ ^s^,, the standard spin correla- 
tion length ^ is dominated by vortices ^ ~ ^i,, while, 
in the low-temperature regime, ^ is dominated by spin- 
waves ^ ~ ^sw Note that, because of the finiteness of 
£,sw, the full spin correlation length is kept finite even 
at T < Ty. The relative importance of vortices and 
spinwaves is interchanged around a crossover tempera- 
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Fig. 2. (Color online) Schematic picture of the decoupling be- 
tween vortices and spinwaves, each characterized by the associ- 
ated correlation length, and ^sw ■ At the Z2-vortex transition 
temperature Ty , the vortex correlation-length diverges, while 
the spinwave correlation length Jstn is kept finite. The behavior 
of the full spin correlation length ^ is also shown. In the low- 
temperature spin-gel state, the spin correlation is dominated by 
spinwave excitations, whereas, in the high-temperature phase, 
the spin correlation is dominated by vortex excitations. There 
exists a crossover temperature Tx (> Ty) at which the relative 
importance of vortices and spinwaves are interchanged. Ty and 
Tx were estimated by the recent Monte Carlo simulation to be 
Ty = 0.285 ± 0.005 and Tx ~ 0.294 (ref. 6). 



turc Tx{> Ty) at which ^„ ~ (Fig.2). Such a de- 
couphng behaviour is in sharp contrast to the behavior 
in other vortex-driven transitions, i.e., the Kosterhtz- 
Thouless (KT) transition of the 2D XY model, where 
there exists only one length scale and the spin correla- 
tion length stays infinite throughout the low temperature 
phase. 

According to a recent theoretical analysis, the di- 
verging behavior of S,v is given by 

A 



£,v exp 



(1) 



where a < 1. The usual thermodynamic macroscopic ob- 
servables such as the specific heat or the magnetic sus- 
ceptibility show only a week essential singularity at T„, 
which makes the direct experimental observation of the 
Z2-vortex transition rather difficult. So far, only an in- 
direct experimental support was reported.^' 

In this paper, we concentrate on the dynamical prop- 
erties of the model in search for the signature of a Z2 
vortex. Our numerical result indicates that the signature 
of a Z2-vortex and a Z2-vortex driven topological transi- 
tion can be captured from the information of dynamical 
spin correlations, specifically from that of the dynamical 
spin structure factor, which can be measured by means 
of inelastic neutron scattering. 

The rest of the paper is organized as follows. In Sec. 2, 
we describe the model and details of our spin-dynamics 
simulation. In Sec. 3, we present the result of our nu- 
merical simulation and analytical calculation, and dis- 
cuss their relation to the Z2-vortex excitation. In Sec. 4, 



we conclude with a summary and implications to exper- 
iments. 

2. Model and Method 

The model we consider is the classical Heisenberg AF 
on the two-dimensional triangular lattice, whose Hamil- 
tonian is given by 



n 



(id) 



(2) 



where J > 0, and the sum is taken over all nearest- 
neighbor pairs. We perform here a spin-dynamics sim- 
ulation of the model. The lattice is a L x L triangular 
lattice of a rhomboidal shape with periodic boundary 
conditions. The dynamics is assumed to obey the classi- 
cal analogue of the Bloch equation. 



_dSi 
I—— 
dt 




X Si, 



(3) 



where the sum is taken over all nearest neighbours. Note 
that this dynamics conserves the total energy as well as 
the total uniform magnetization M — Si . The tem- 
perature effect is taken into account via initial spin con- 
figurations, which are generated from the standard equi- 
librium Monte Carlo (MC) simulation at the temperature 
T. Dynamical observables are calculated by averaging 
over 500 — 1000 independent spin-dynamics simulations 
starting from different initial configurations prepared by 
equilibrium MC runs. 

In order to prepare initial spin configurations for our 
spin-dynamics simulations, we generate equilibrium spin 
configurations at a given temperature T by first perform- 
ing a hybrid Monte Carlo simulation which consists of 
the heat-bath updating and the over-relaxation updat- 
ing. One hybrid MC step (MCS) consists of one heat-bath 
and ten over-relaxation sweeps. At each temperature, we 
generate such equilibrium spin configurations, which are 
taken from our MC run every 10'^ MCS after discarding 
10^ MCS for equilibration. 

In our spin-dynamics simulation, we employ a recently 
developed second-order decomposition method to solve 
equation (3) numerically.^''' We set a time step 5t = 
Q.OlhJ^^, and the integration of equation (3) is carried 
out to tmax, typically of the order of tmax — 800hJ~^. 
We checked the accuracy of our simulation by compar- 
ing the obtained data at several temperatures with those 
obtained with a finer time step of St = 0.001?iJ~^, and 
with those obtained by using the fourth-order Runge- 
Kutta method with a time step of St = O.OlftJ"^. 

The dynamical spin structure factor is defined by 

S{q,uj) = J dr J dt (^So{0) ■ Srit)"^ cxp[-i{ujt + q ■ r)]. 

(4) 

where uj is the angular frequency, q is the wavevector, 
and (• • ■ ) denotes a thermal average. 

The dynamical spin structure factor is calculated via 
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the formula, 



S{q, 



(5) 



where (• • • ) represents an average over initial spin config- 
urations, and Sq{uj) is the space-time Fourier transform 
of the spin variable given by 

Sg{Lj) = j dtY^ Si{t) exp[-i(q • n + cot)]. (6) 

i 

In order to decrease statistical errors, the dynamical spin 
structure factor is symmetrized as 



S{q.i 



(7) 



The central component of S{q,w) is often subject 
to pronounced finite-size effect. We have checked the 
system-size dependence of S{q,uj) by comparing those 
of L = 192, 384, 768 at each temperature for the triangu- 
lar case, to confirm that the data shown here are close 
enough to the bulk ones. 

3. Results 

In Fig.Sa, we show the frequency dependence of the 
dynamical spin structure factor S{q, uj) at a wavevector 
q = (47r/3-H47r/192, 0), slightly off the K-point. Note that 
the perfect 120-degrees structure yields the magnetic 
Bragg peak at the K-point, qx = (47r/3,0). The tem- 
perature range studied is below around the specific-heat 
peak temperature Tpeak — 0.32, including both above 
and below Ty ~ 0.285. As can be seen from the figure, 
at a temperature T = 0.280 slightly below Ty, S{q,u)) 
exhibits only two side peaks corresponding to spinwave 
excitations. On increasing the temperature, in addition 
to the spinwave peaks, a central peak centered at w = 
appears and grows. The central peak and the spinwave 
peaks coexist in the range 0.290 < T" < 0.305, whereas, 
at higher temperatures, these peaks are merged into a 
single broad peak at w = with spinwave shoulders re- 
maining at finite ui. At T = 0.305, the intensity of the 
central peak is about 1% of that at the K point. The 
width of the central peak is of the same amount as that 
of the spinwave peak, showing no pronounced tendency 
to sharpen toward T — >■ T„. Note that the onset of the 
central component in S(q,uj) is fully correlated to the 
vortex transition at T = T„ ~ 0.285. This correspon- 
dence suggests that the observed central peak might be 
originated from the contribution of free Z2 vortices. 

In order to further clarify the situation, we also com- 
pute S{q, uj) for the unfrustrated 2D Heisenberg AF on a 
square lattice, a model which does not bear a Z2 vortex. 
The result is shown in Fig. 3b for the temperatures below 
around the specific-heat peak temperature Tpeak — 0.68. 
In sharp contrast to the triangiilar case, no central peak 
is observed in the low-temperature regime where spin- 
wave peaks exist. This absence of the central peak co- 
existing with the spinwave peaks supports the view that 
the central peak observed in the triangular case is due to 
free Z2 vortices. 

We next perform a simplified theoretical analysis of the 
dynamical spin structure factor based on the spinwave- 



vortex decoupling picture. We start with a simple as- 
sumption of the factorization of the dynamical spin-spin 
correlation function C{r,t) into the vortex part and the 
spinwave part 

C{r,t) = (5o(0) • Sr{t)) ~ Cy{r,t)CsUr,t), (8) 

where Cs^(r,t) is a dynamical spin correlation func- 
tion due to the spinwave excitations and Cv{r,t) is the 
one due to the Z2-vortex excitations. We assume that 

the space Fourier transform of Csjj,>(T",t) is governed by 
the spinwcLve mode with a dispersion uJsw 

(q) and with a 

damping Tsw{q) as 

Csw{q,t) = Ssn,{q)exp[-t/Tsu,{q) ± iuJsw{q)t\, (9) 

where Sswiq) is the static spin structure factor due 
to spinwaves. We also assume that the corresponding 
vortex-part is given by a simple exponential form as 



Cv{q,t) = Sv{q)exp[-t/Tv], 



(10) 



where Sv{q) is the static spin structure factor due to 
vortices, and t„ is assumed to be q-independent. From 
eqs. (8)-(10), S{q,co) is obtained as the convolution of 
the spinwave part and the vortex part as 



S{q 



,w) = j dq'- 



2Tiq - q')SUq - q')Sv{q') 
+ [T{q-q')nuj±uJsy,{q-q')Y 



(11) 

where the reduced correlation time is given by l/r^q) = 
i/Tsw{q) + l/^L.. The peak of Ssw{q) is located at the 
K point, while Sv{q) is expected to be peaked at g = 
0. Thus, S{q,uj) typically exhibit two peaks, one at 
oj = ujsw{q) and the other at w = 0Jsw{qK)- Since 
^sw{,qK) = 0, the latter corresponds to the central peak. 
At the Z2-vortex transition T = Ty, the vortex correla- 
tion length ^„ diverges, and in the low temperature spin- 
gel state, Sy{q) oc S{q). Thus, in the spin-gel state below 
Ty, S{q,u}) is determined solely by spinwaves without a 
central-peak structure. 

In Fig. 3c, we show a typical form of S{q, lu) above Ty 
calculated from eq. (11) under several assumptions given 
in the caption of Fig. 3. As can be seen from the figure, 
the observed shape of the dynamical spin structure factor 
is well reproduced by (11), at least qualitatively, as long 
as T{q) is taken to be moderately small. This observation 
gives further support to the view that the central peak 
of S{q, u) is originated from free Z2 vortices. 

Another explanation of the central peak of S{q, uj) was 
suggested by Mertens et al in terms of the KT vortices 
of the 2D XY model. Based on a dilute vortex-gas 
picture, these authors estimated the free-vortex contri- 
bution to S{q, w) in the form of squared Lorenzian.^^' 
The data of several spin-dynamics simulations on the 2D 
XY model were interpreted in terms of this theory. ^^"^-'^^ 
In our observation, however, the motion of Z2 vortices 
are incompatible with such a dilute gas picture. The Z2 
vortices of the present model behave diffusively so that 
the root-mean-squared displacement of the free vortices 
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Fig. 3. (Color online) (a) The frequency dependence of the dynamical structure factor of the triangular-lattice Heisenberg antiferro- 
magnet at a wavevecotr q = (47r/3 + 47r/192,0) near the K point for the temperatures below around the specific-heat maximum 
temperature T < Tp^^k — 0.32 including both below and above T„ ~ 0.285. The lattice size is L = 768. For better visualization, 
baselines of the data at different temperatures are shifted vertically (dotted lines), (b) The frequency dependence of the dynamical 
structure factor of the square-lattice Heisenberg antiferromagnet at a wavevecotr q = {n — 27r/128,7r — 27r/128) for the temperatures 
below around the specific-heat maximum temperature T < T^^^k — 0.68. The lattice size is L = 256. (c) A typical pattern of the 
dynamical structure factor calculated from eq. (11) by assuming Lorenzian forms for Ssw{q) and Sv{q), Sswiq) oc 1/[1 + (q — qic)^] 
and Sv{q) oc 1/[1 + ^v1^]t with a linear dispersion relation, LUs-w{q) = Vsw\q — qK\ and neglecting the q-dependence of 'r(q). The 
parameters in cq. (11) are = 250, (,sw = 150, q = (47r/3 -|- 47r/192, 0), T(q) = 40 and Vaw = 1.36. The central peak arises when t is 
moderately small, l/uJsw{q) ^ t < S^sw/vsw 




hw/J hoj/J 



Fig. 4. (Color online) (a) The frequency dependence of the dynamical spin structure factor of the triangular-lattice Heisenberg antifer- 
romagnet at a wavevecotr q = (0, 27r/\/3) corresponding to the M point for the temperatures as in Fig. 3a. The lattice size is L = 768. 
For better visualization, baselines of the data at different temperatures are shifted vertically (dotted lines). As a guide to the eye, the 
central component referred in the main text is indicated by the dashed lines, (b) The frequency dependence of the dynamical structure 
factor of the square-lattice Heisenberg antiferromagnet at a wavevecotr q = (0, tt) for the temperatures as in Fig. 3b. The lattice size 
is L = 256. 
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is proportional to the square root of the elapsed time t 
as -y/ {dr'^{t)) oc ^/t. Hence, the origin of the central peak 
observed in the present model probably differs from the 
one proposed by Mertens et al for the KT vortices. 

The central peak of S{q, lo) is not limited to the vicin- 
ity of the K point. A broad central component also ap- 
pears even at wavevectors away from the K point, though 
its intensity is significantly reduced and its width is 
much broadened. In Fig. 4a, we show the oj-dependence 
of S{q,ijo) at the M point, q = Qm = (0, 27r/A/3). In ad- 
dition to the spinwave side peaks, a broad but distinct 
central component centered at w = exists. Note that 
the width of the central component is larger than that 
near the K point by an order of magnitude, while its in- 
tensity is weaker than that near the K point by a factor 
of 10^^. The intensity of the central component tends to 
be decreased with decreasing the temperature as in the 
vicinity of the K point. 

For comparison, we show in Fig. 4b S{q, lj) of the 
square-lattice Heisenberg AF at the zone boundary, q = 
(0, tt). The computed S{q,uj) hardly exhibits a central 
peak. Though a faint central component is discernible, 
its intensity is significantly weaker than that of the spin- 
wave peak, and the difference from the triangular case 
shown in Fig. 4a is distinct. 

We then ascribe the origin of the broad central compo- 
nent observed away from the K point to the contribution 
of the Z2 vortex-pair excitations. Near the K point, the 
central peak of S{q, uj) predominantly probes the dynam- 
ics of free Z2 vortices, while, as one moves away from the 
K point, the central component tends to probe the dy- 
namics of more tightly bound Z2-vortex pairs. Our obser- 
vation that the width of the central component becomes 
significantly larger when one moves away from the K 
point (see Fig. 3a and Fig. 4a) means that the time scale of 
underlying dynamics becomes considerably faster. This 
seems fully consistent with our identification above, since 
tightly bound pairs tend to yield faster dynamics asso- 
ciated with pair annihilation and creation than that of 
free vortices. Another feature to be noticed is that, at 
the q points away from the K point, e.g., at the M point, 
the central component of S{q,Lu) tends to persist even 
at T < Ty. This observation also supports our interpre- 
tation above, since the vortex pair can stably exist even 
below Ty unlike the free vortex. 

In order to illustrate the broad central component ob- 
served for a wide range of wavevecotrs, we show in Fig. 5 a 
the intensity map of 5'(<7, uj) in the entire q-planc at 
huj/J ~ 1 at a temperature slightly above Ty, T — 0.31. 
The effect of the central components is observed as the 
intensities near the Brillouin zone boundaries, which ap- 
pear in addition to bright ring-like intensities surround- 
ing the K points corresponding to spinwave excitations. 
We note that such intensities near the zone boundaries 
cannot be seen in the corresponding intensity map of the 
square-lattice Heisenberg AF shown in Fig. 5b. 

4. Summary and Discussion 

From all these observations, we conclude that the cen- 
tral peak of the dynamical spin structure factor S{q, uj) 
is originated from the Z2-vortex excitation. The bases 



(a) 

Triangular lattice {T=0.31, = l.OJ) 




(b) 



Square lattice (T=0.66, fvjj = l.OJ) 




Fig. 5. (Color online) (a),(b) the intensity maps of the dynamical 
spin structure factor in the wavevecotr space with hu) = J, for the 
case of the triangular-lattice Heisenberg antiferromagnet with 
L = 192 at a temperature T = 0.31 (a), and for the case of 
the square-lattice Heisenberg antiferromagnet with L = 256 at a 
temperature T = 0.66 (b). 

of our conclusion can be summarized as: i) Near the K 
point, the central peak appears only &i T > Ty, increas- 
ing its intensity with increasing the temperature; ii) the 
central peak does not appear in the unfrustrated square- 
lattice Heisenberg AF, a model which does not sustain 
the Z2 vortex; and iii) the central-peak structure similar 
to the one observed in our simulations can be reproduced 
by a simplified theoretical model based on the spinwave- 
vortex decoupling picture. 

In contrast to the information of the standard spin 
correlation length ^, which is obtainable via the elas- 
tic neutron scattering measurements, the information of 
the vortex correlation length £^y, that is the length scale 
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diverging at the vortex-transition point T = T^, is hard 
to obtain. In principle, however, appropriate information 
about the central peak enables one to deduce ■ On the 
basis of (11) together with the Lorenzian forms of Ssw{q) 
and Sv{q), the ratio of the central-peak intensity Ic{q) to 
the spinwave one Isw {q) near the K point At q ^ qk + Sq 
is approximately obtained as Id Isw ^ S,sw/S,v By com- 
bining it with the information of the standard spin cor- 
relation length ^ ^ S.swS.v/iS.sw we can estimate £,y 
together with 

Experimentally, the dynamical spin structure factor 
can be measured via inelastic neutron scattering. As 
mentioned, candidate triangular-lattice AFs possessing 
a Z2 vortex excitation might be NaCr02, NiGa2S4 and 
organic compounds like k-(BEDT-TTF)2Cu2(CN)3 or 
EtMe3Sb[Pd(dmit)2]2- Unfortunately, the organic com- 
pounds are not suited to neutrons so that neutron- 
scattering measurements seem unrealistic for these ma- 
terials at the present stage. Then, promising candidate 
materials for the Z2-vortex detection might be NaCr02 
and NiGa2S4. These materials exhibit a transition- 
like anomaly at T ~ 30K (NaCr02) and T ~ 8.5K 
(NiGa2S4), slightly below the specific- heat peak tem- 
perature. Often, only powder samples are available for 
these materials. Under such circumstances, we also com- 
pute the angle-averaged dynamical spin structure factor 
S{\q\,uj) to examine how S{\q\,uj) looks like for powder 
samples, and the result is shown in Fig. 6a and 6b for 
both cases of the |q|-value close to |g| ~ \qK\ (a), and 
\q\ ~ \qM\ (b). In both cases, the central peak is clearly 
observable together with the spinwave peaks. Thus, in- 
elastic neutron scattering experiments for powder sam- 
ples would be sufficient for observing the signature of a 
Z2 vortex and a Z2-vortex driven topological transition. 
This is a nice aspect for practical purpose, since one can 
expect the intensity gain for powder measurements. 

Inelastic neutron-scattering measurements as sug- 
gested here might provide a useful information in under- 
standing the spin-liquid-like behavior and the transition- 
like anomaly observed in triangular-lattice Hcisenbcrg 
AFs NaCr02 and NiGa2S4, and the nature of novel ex- 
citations in geometrically frustrated magnets in general. 
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Fig. 6. (Color online) The frequency dependence of the angle-averaged dynamical spin structure factor of the triangular-lattice Heisen- 
berg antiferromagnet for the temperatures as in Fig. 3a and Fig. 4a. The data should correspond to the neutron-scattering pattern 
obtained from powder samples. Panels (a) and (b) correspond to wavevectors satisfying \qK\ + 47r/192 < \q\ < IqkI +67r/192 (a), and 
IqmI - 277/384 < |<7| < |<7mI + 27r/384 (b). 



